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We propose to use the ratio of the pion-proton correlation 
functions evaluated under different conditions to study the 
relative space-time asymmetries in pion and proton emission 
(pion and nucleon source relative shifts) in high energy heavy 
ion collision. We address the question of the non-central col- 
lisions, where the sources can be shifted spatially both in the 
longitudinal and in the transverse directions in the reaction 
plane. We use the RQMD event generator to illustrate the 
effect and the technique. 

PACS number: 25.75. Gz 



The importance of the study of the space-time struc- 
ture of particle emission in the heavy ion collisions has 
been emphasized repeatedly. In general the question 
of the corresponding measurements splits into two: the 
measurements of the extensions (sizes) of the emission 
zones of each particular particle species and the mea- 
surement of the relative location of the sources of dif- 
ferent particles. The measurements of the size of the 
effective source could be addressed in principle by the 
two (identical) particle interferometry. At present the 
particle sources are studied intensively by this method. 
On the other hand, the question of the relative space and 
time asymmetry in the production of different particles 
(source shifts) remains almost unexplored, although this 
question is very important for an interpretation of the 
experimental data within different models (for example, 
the ones which require particle thermalization, chemical 
equilibration, etc.). 

The important first step in the investigation of the 
source shifts in space and time was done in Ref . [Q , where 
the non-identical particle correlation functions were pro- 
posed as a tool for a study of time delays in emission 
of different particle. This observation was based on the 
analytical calculations of the final state interaction 
contribution to the correlation function. Recently, more 
detailed investigation of the possibilities provided by the 
correlations of two non-identical particles for the study of 
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asymmetries in particle production was attempted in [Q . 

Here, we use the n^p correlations to study the asym- 
metries in the pion and nucleon (proton) production in 
non-central nucleus- nucleus collisions at high energies. 
To be specific, we investigate the particle production 
in the rapidity region close to the projectile rapidity in 
Au-|-Au collisions at the AGS energies. We select the 
rapidity region where proton directed flow is most pro- 
nounced 1^ , and where it is natural to expect that proton 
and pion sources could be shifted relative to each other 
not only in the longitudinal (z) direction but also in the 
flow direction (along or opposite to the impact parameter 
vector, below, +x or —x directions). 

The two-pion interferometry of the non-central colli- 
sion and its relation to anisotropic flow was addressed 
earlier in |c|. In particular, it was noted that the pion 
source looks different from the different directions with 
respect to the reaction plane angle, in part due to the 
fact that the pion source is screened by nucleons from 
the direction of nucleon flow. It implies that pion source 
is effectively shifted with respect to the proton source. 
This observation was supported by the analysis of the 
pion source shape in the RQMD generated events. The 
recent experimental results on anisotropic flow of pions 
and nucleons Q confirm the picture. The observed flow 
of protons and pions looks very different, and this dif- 
ference agrees with the assumption of pion screening by 
nucleons. What is even more important that pions of dif- 
ferent charge at very low transverse momenta flow in the 
opposite directions. This fact finds natural explanations 
in Coulomb interaction with the protons, assuming that 
the pion and proton sources are shifted in space. In the 
current paper we show how such shifts could be detected 
and measured experimentally using pion-proton correla- 
tions, determined mainly by the two-particle final state 
interaction. 

Referring for the quantitative description of the ques- 
tion of the final state interaction contribution to the cor- 
relation function to the original papers here we 
discuss it qualitatively. Such a discussion will help us to 
explain/justify the techniques used below. 

Any two-particle (both identical and non-identical) 
correlation function is mainly sensitive to the distance 
between particles at the moment when the second parti- 
cle is produced. Such a separation can be written as 
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fi2 « (ri - r2) - V{ti - t2), 



(1) 



where we introduced the notations r, t for the space 
and time position of particle production, and V — 
(Pi + P2)/(-Ei + E2) for the pair velocity. The differ- 
ence between identical and non-identical particle corre- 
lations (except the quantum statistics effects, which are 
not essential for the current study) comes from the fact 
that the correlation function in the case of non-identical 
particles in general depends on the orientation of the rel- 
ative velocity v = pi/i?i — P2/-E2 (that is, the correlation 
function is different for the cases of v and — v), while for 
identical particles such a dependence is washed out due 
to the symmetry of the system. 

The dependence of the non-identical two-particle cor- 
relation function on the relative velocity becomes clear 
if one recollects that the non-identical two-particle cor- 
relations are mostly due to the two body interaction in 
the final state. The impact of the two body final state 
interaction obviously is different for particles (at the in- 
stant of both of them are created) moving, in the particle 
pair rest frame, toward each other or moving in the op- 
posite directions. If the centers of the effective sources 
of two particles are shifted, then the dependence of the 
correlation function on the orientation of the relative ve- 
locity with respect to the shift becomes obvious. Suppose 
that the particles of type "2" are produced on average at 
Z2 > zi (for simplicity assume also that ^2 ~ ^i)- Then, 
selecting the pairs with vi^z > V2,z one would select the 
pairs which move on average toward each other (in the 
pair rest frame) and the corresponding final state interac- 
tion would be stronger. (In this simplified consideration 
we also neglect the possible differences between particle 
velocities at the moment of the production and the mea- 
sured ones). Depending on the type of final state interac- 
tion the correlation function would decrease or increase. 
Studying the orientation of the relative velocity to which 
the correlation function is the most sensitive one could 
find the orientation of the shift in space. 

It is natural to study the two non-identical particle 
(tt^p) correlations in the pair rest frame and as a func- 
tion of k*, the (half of) particle relative momentum in 
this system (k* = Ptt = — Pp)- In the region of small val- 
ues of k* (much less than the inverse size of the particle 
source) the final state interaction of two charged parti- 
cle is dominated by Coulomb interaction. In this case 
the correlation function can be written in the analytical 
form l^-Ql, which makes it possible to estimate the ra- 
tio of the correlation functions evaluated under different 
conditions. If we denote the correlation functions for two 
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different cases, for example, k* > and k* < 0, as R 
and i?,- \ then in the limit of small fc*, 
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a + 2(r*) + 2(r*k*/fc*)(+) 
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where a is the Bohr radius (for the Tr^p system a « 
±222 fm), r* is the relative separation in particle pro- 
duction points. Then, (r*) is the shift between sources 
in the pair rest frame, the quantity of interest. It should 
be mentioned, that in the derivation of Eq. (|^) it is as- 
sumed that (r*) ^ |a|, and that |Re/| <C |a|, where / is 
the strong s-wave pion-proton scattering amplitude. The 
average in formula (^) is taken over all pairs from the 
corresponding subsample. It was also used that in the 
limit of k* < (pt), (r*)(k*/fc*)(±) = ±(i-*>^/2 for the 
above defined cuts on k* . 

For the current study we use the RQMD vl.08 § 
event generator to simulate Au-|-Au collision at 
11.4 GeV/nucleon. We select particles in the rapidity 
region 2.8 < yiab < 3.2 and within relatively narrow sec- 
tor in the azimuthal space, always such that Px > ^ and 
\Py/Px\ < 0.5. We create two event subsamples in accor- 
dance to the orientation of the reaction plane, "'5,. = 0" 
(in this case nucleons flow in the positive x direction) 
and "^f^ = tt" subsamples. The relative shifts between 
pion and proton sources for both cases are presented in 
Table 1 as calculated both in the center of mass system of 
colliding nuclei and in the particle pair rest frame. In the 
RQMD generated events we observe strong momentum- 
position correlations in particle production (especially for 
pions); the values presented in the table correspond to 
the average of the production points over the particles 
which contribute to the pairs with k* < 15 MeV. This 
restriction limits significantly the effective transverse mo- 
mentum range of pions, decreasing its average value in 
the sample. For both subsamples the average compo- 
nents of the pair velocity in the center of mass system of 
colliding nuclei are Vz ~ 0.89 and Vx ~ 0.17. 

We analyze approximately 200K pairs in each of the 
subsamples, the typical statistics for modern experi- 
ments. Our goal is to show the sensitivity of the irp 
correlation function to the shifts presented in Table 1 

Below we consider only Tr~^p correlation function 
(which decreases with the strength of the interaction), 
the results for tt~p case would be very similar taking 
into account the change in sign of Coulomb interaction, 
which dominates at low values of particle relative veloc- 
ity. The correlation functions are calculated in accor- 
dance to Lednicky-Lyuboshitz j^] formulae taking into 
account both Coulomb as well as the strong interaction 
in the final state. In Fig. 1 we show the correlation func- 
tions calculated for the "^^r — 0" event subsample for 
two sets of cuts A;* > (fc* < 0) and k* > {k* < 0). 
Note that these cuts approximately correspond to the 
cuts in the laboratory system > Vp^x ivTT,x < Vp^x) 
and Vt^^z > Vp^z {vtt,z < ^'p,z), respectively. The corre- 
spondence would be strict if one selects the coordinate 
system with one of the axes parallel to the pair velocity. 
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The ratios of the correlation functions are also shown on 
the same figure. One can see that the correlations are 
stronger in the cases of fcj > and k* > 0, which clearly 
indicates that the proton source is shifted relative to the 
pion source to positive z* and positive x* values. The 
magnitude of the difference in the cases of cuts on /c* 
and k* shows that the shift in z direction is larger than 
the one in x direction. 

In Fig. 2 we show the analogous plots for the "'I'r ~ tt" 
subsample. Note the difference in Fig. 1 an Fig. 2 con- 
cerning the cut on k* . It reflects the fact that the relative 
shift between pion and proton sources has changed the 
sign (and its value), and the proton source is now to the 
"left" relative to the pion source. The changes in the 
correlation function ("A:*" cuts) although small in mag- 
nitude are significant statistically, what is seen in Fig. 3, 
where we show the correlation functions at a different 
scale. 

Qualitatively the results of the tt~^p correlation anal- 
ysis show unambiguously the correct relative location of 
pion and proton sources. Quantitatively the agreement 
with the approximate formula (^) is also very good. In 
accordance to (^), each one fermi of the shift leads to 
the change of approximately 0.9% in the ratio of the cor- 
responding correlation functions at small values of k*. 
If one compares the numbers from the Table I with the 
values of the correlation function ratios extrapolated to 
k* = from Figs. 1-3, a good quantitative agreement will 
be found. From the Table I one would expect the follow- 
ing values for the ratios of the correlation functions at 
small values of k*: in the case of "^f^ — 0" the ratio 
Ri~^/Ri'^^ ~ 1.11 and R^^^ /R^"^^ « 1.05; in the case 
of = tt", Ri'^/Ri'^^ « 1.09 and i?l"V-R^^^ ~ 0-99. 
The correlation function ratios shown in Fig. 1-3 agree 
with these numbers quite well. 

In the current paper we have analyzed the correlation 
function in the particle pair rest frame. The source shifts 
which were extracted are the ones in this frame. Al- 
though, in principle, the correlation function is sensitive 
to the time shift between the sources in this system, 
the sensitivity is very weak and can be neglected under 
the condition of ^ fir*^ (fi is the reduced mass). 
Physically, one of the reasons for this is the very small 
particle relative velocity in this frame v* « k*//z; then 
the relative space separation between particles changes 
on average very little during the time At ~ t* provided 
that k*/fi\At\ < |{r*)|. 

If one considers the space and time shifts in the labo- 
ratory system (or in the center of mass of colliding nu- 
clei system) then the correlation function depends on all 
four shifts, but all four of them cannot be extracted from 
3 possibly measured ratios of the correlation functions. 
It is exactly the same problem as observed in the two- 
particle interferometry. To extract all shifts in this case 
one needs to make additional assumptions, for example, 



that the shifts do not depend on the velocity of the pair. 
Such questions are beyond the scope of the current pub- 
lications. 

We have shown that the important information on the 
relative shifts between pion and proton sources can be 
obtained by comparing the wp correlation functions eval- 
uated at different conditions. The good quantitative re- 
sults achieved in the application of the method to the 
RQMD generated events give a hope that the same tech- 
nique can be successfully applied to the data. 
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TABLE 1. The mean values of spatial and temporal shifts 
(in fm) of pion and proton sources for two different orientations 
of the reaction plane in the system of the center of mass of collid- 
ing nuclei (upper half) and in the pair rest frame (lower half). 
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FIG. 1. The correlation functions -R^ti (k*c,z > 0) and 
^x,l {kx,z < 0) for the event subsample "vPr = 0". 
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FIG. 2. The same as Fig. 1 for the event subsample 
=7r". 
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FIG. 3. The ratios of the correlation function ^ and 
for two event subsamples "*r = 0" and "*r = tt" . 
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